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ABSTRACT 

We present a search for galaxies at 7.6 < z < 9.8 using the latest HST WFC3 near- 
infrared data, based on the Lyman-break technique. We search for galaxies which have 
large (Y—J) colours (the "K-drops") on account of the Lyman-a forest absorption, and 
with (J — H) colours inconsistent with being low-redshift contaminants. We identify 
24 candidates at redshift z « 8 — 9 (15 are robust and a further 9 more marginal but 
consistent with being high redshift) over an area of w 50 square arcminutes. Previous 
searches for F-drops with WFC3 have focussed only on the Hubble Ultra Deep Field 
(HUDF), and our larger survey (involving two other nearby deep fields and a wider 
area survey) has trebelled the number of robust F-drop candidates. For the first time, 
we have sufficient z w 8 — 9 galaxies to fit both <j>* and M* of the UV Schechter 
luminosity function. There is evidence for evolution in this luminosity function from 
z = 6 — 7 to z = 8 — 9, in the sense that there are fewer UV-bright galaxies at z w 8 — 9, 
consistent with an evolution mainly in M* . The candidate z s=s 8 — 9 galaxies we detect 
have insufficient ionizing flux to reionize the Universe, and it is probable that galaxies 
below our detection limit provide a significant UV contribution. The faint-end slope, a, 
is not well constrained. However, adopting a similiar faint-end slope to that determined 
at z = 3 — 6 (a = — 1.7) and a Salpeter initial mass function, then the ionizing photon 
budget still falls short if / esc < 0.5, even integrating down to Mjjv = —8. A steeper 
faint end slope or a low-metallicity population (or a top-heavy IMF) might still provide 
sufficient photons for star-forming galaxies to reionize the Universe, but confirmation 
of this might have to await the James Webb Space Telescope. 

Key words: galaxies: evolution galaxies: formation galaxies: starburst galaxies: 
high- redshift ultraviolet: galaxies 



1 INTRODUCTION 

The z « 8 epoch is cosmologically very interesting: the 
Gunn-Peterson effect (Gunn & Peterson 1965, Scheuer 
1965), the near total absorption of the continuum flux at 
wavelengths shorter than that of the Lyman-a line due to 
a significant neutral hydrogen fraction in the inter galactic 
medium (IGM), has been observed at z > 6.3 in SDSS QSO 
spectra (Becker et al. 2001; Fan et al. 2001, 2006). This sug- 
gests that z w 6 lies at the end of the Epoch of Reionization, 
whose mid-point may have occurred at z « 11, according to 
latest results from WMAP (Dunkley et al. 2009). 

An outstanding problem, however, is what sources were 
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responsible for the reionization of the Universe, and when 
exactly this occurred. There is evidence for old stellar pop- 
ulations in some z ~ 4 — 6 galaxies from Balmer break mea- 
surements in Spitzer/IRAC imaging (Eyles et al. 2005, 2007; 
Stark et al. 2007, 2009), implying star formation commenced 
at even earlier times. Hence it is reasonable to consider the 
UV photons from this star formation as a possible cause 
of reionization. Age and stellar mass determinations of the 
above mentioned stellar populations are affected by many 
uncertainties, so is important to directly look for star for- 
mation at z > 7 to determine whether star-forming galaxies 
at these epochs can indeed provide the Lyman continuum 
photons required for reionization. 

In recent years, observations of high redshift universe 
(z > 6) have become possible. Deep imaging surveys with 
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the Hubble Space Telescope (HST) and large ground based 
telescopes have made the discovery of z ~ 6 galaxies al- 
most routine. Some of those searches (Bunker et al. 2004, 
Yan & Windhorst 2004, Stanway, Bunker & McMahon 2003, 
Bouwens et al. 2006, 2007, Oesch et al. 2007, Yoshida et al. 
2006) rely on the Lyman break galaxy (LBG) technique, 
initially used by Steidel and collaborators (Steidel et al. 
1996) to identify galaxies at z « 3 through the large ab- 
sorption produced by the intervening Lyman-a forest clouds 
and the Lyman limit. Until recently, working at higher red- 
shift (z > 7) was challenging. Studies were limited to very 
small deep fields observed from space (e.g. Bouwens et al. 
2008), or to extremely shallow wide-area surveys from the 
ground (e.g., Stanway et al. 2008, Hickcy et al. 2010). The 
new Wide Field Camera 3 (WFC3) instrument, installed on 
HST in May 2009, allowed this technique to be more ef- 
fectively applied to z « 7 — 10, thanks to its near-infrared 
channel with significantly larger field and better sensitivity 
than the previous-generation NICMOS instrument. Using 
WFC3 broad-band filters at 1.0, 1.25 and 1.6/im (the Y-, 
J- and //-bands) and targetting fields with existing deep 
Advanced Camera for Surveys (ACS) data it is possible to 
identify optical "drop-outs", objects seen only in the WFC3 
infrared images but not in the optical ones. These are can- 
didate z > 7 galaxies. Searching for objects with no flux at 
1.0/im and below (Y-band drop-outs, or "Y-drops") could 
lead to the discovery of z « 8 galaxies. Deep optical images 
of the observed fields are still necessary to "clean up" the list 
of candidates, because, as we will see later, optical detections 
are useful in ruling out many lower-redshift contaminants. 

The past few months have seen several papers present- 
ing high redshift galaxy candidates from HST/WFC3 imag- 
ing of the Hubble Ultra Deep Field (HUDF; Bunker et al. 
2010, McLure et al. 2010, Oesch et al. 2010, Bouwens et al. 
2010a, Yan et al. 2010, Finkelstein et al. 2010). In the HUDF, 
~ 10 z'-drops (z w 7) have been found, along with ~ 5 Y- 
drops (z « 8). Spectroscopic confirmation of these candi- 
dates in the HUDF will be extremely challenging, as they 
have magnitudes J > 26.5 (for the z'-drops) and J > 28.0 
(for the Y-drops). What is needed are larger samples over 
wider areas, which might yield rarer but brighter candidates 
more suitable for spectroscopic follow-up. Such follow-up is 
important to test the validity of the Y-drop selection tech- 
nique, and address the contaminant fraction, as well as ex- 
ploring the physics of star forming galaxies at z > 8 (in par- 
ticular whether Lyman-a emerges during the Gunn-Peterson 
absorption era). Increasing the survey area of the WFC3 
LBG searches will also improve the statistics (and hence 
the rest-UV luminosity function constraints), and we have 
started to do this by searching for z'-drops in the larger- 
area Early Release Science WFC3 images of some of the 
GOODS-South field (Wilkins et al. 2010a) and expanding 
this to include two other deep flanking fields (UDF-P12 & 
UDF-P34) close to the HUDF (Wilkins et al. 2010b), which 
has increased the number of robust z'-drops from ~ 10 to 
~ 40. In this paper we use our new reductions of the ERS 
GOODS-South, and UDF-P12 & UDF-P34 to search for Y- 
drops at z w 8 — 9. In Bunker et al. (2010) we presented our 
preliminary list of Y-drops in the HUDF, and here we also 
re-analyse this field using a more recent data reduction. 

The evolution of the rest-UV luminosity function is key 
to both understanding the star formation history of the Uni- 



verse, and also to address the role of star-forming galaxies 
in reionization. There seems to be strong evolution in the 
UV luminosity function (UVLF) up to redshift up to z « 6 
(e.g. Stanway, Bunker & McMahon 2003), and recent studies 
(Bunker et al. 2010; Wilkins et al. 2010b; Oesch et al. 2010) 
seem to show that this evolution continues up to z « 7, al- 
though based on small-number statistics. Our goal is to push 
the measurement of the UVLF further back in cosmic time 
by assembling a statistically-significant sample of probable 
z « 8 — 9 galaxies. From this we can address the evolution 
of the star formation rate density, and the ionizing photon 
budget. 

This paper is organised as follows: in Section 2 we out- 
line the HST observations with WFC3 and the data reduc- 
tion, and in Section 3 we describe our colour selection to 
recover high-redshift Lyman break galaxies, and compare 
our sample with those from other studies. In Section 4 with 
discuss the evolution of the star formation rate density, and 
the implications for reionization, derived from the luminos- 
ity function we infer at z ~ 8. Our conclusions are pre- 
sented in Section 5. Throughout, we adopt the standard 
concordance cosmology of Qm = 0.3, = 0.7 and use 
Ho = 70kms _ Mpc -1 . All magnitudes are on the AB sys- 
tem (Okc & Gunn 1983). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

In this paper we analyse images from WFC3 on HST taken in 
the near-infrared Y-, J- and //-bands. The data come from 
two different HST programs, both covering areas within the 
GOODS-South field (Giavalisco et al. 2004). The HST Trea- 
sury programme GO- 11 563 (P.I. G. Illingworth) covers the 
HUDF and two nearby deep flanking fields (UDF-P12 and 
UDF-P34, also referred to as HUDF05-01 and HUDF05-02 
in programme GO-11563). These flanking fields were imaged 
by the Advance Camera for Surveys (ACS) on HST in v-, 
i'- and z'-bands during 2005-6 in parallel with deep HST 
NICMOS NIC3 observations of the original UDF as part 
of program GO-10632 (P.I. M. Stiavelli). In Bunker et al. 
(2010) we analysed the single- WFC3-pointing HUDF data 
obtained soon-after the commissioning of WFC3, and in this 
paper we study the two new deep WFC3 pointings on the 
two deep flanking fields, and reanalyse the UDF data using 
more recent on-orbit calibration of the detector. Addition- 
ally, we analyse the Early Release Science (ERS) program 
GO/DD-11359 (P.I. R. O'Connell) data, covering ten over- 
lapping pointings with two orbits in each filter. An analysis 
of the first 6 pointings for z-drops at z « 7 was presented in 
Wilkins et al. (2010a), with the full ERS mosaic and UDF- 
P12&P34 flanking fields used to select z-drops in Wilkins 
et al. (2010b). 

The infrared channel of WFC3 was used, which is a 
Teledyne 1014 x 1014 pixel HgCdTe detector (a 10-pixel 
strip on the edge is not illuminated by sky and used for 
pedestal estimation), with a field of view of 123" x 136". 
Filters used in the two programs are the same for J and H 
bands (F125W and F160W), while the ERS images use a Y- 
band filter (F098M) which covers only the blue side of the 
wider F105W filter used in the UDF and flanking field im- 
ages. The data were taken in "MULTIACCUM" mode using 
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WFC3 exposure times, in ksec (number of exposures). 



Field ID 


F-band a 


J-band 


iJ-band 


J 6a 


J 7cr 


HUDF 


39.3 (28) 


44.9 (32) 


78.6 (56) 


28.51 


28.34 


P34 


28.1 (20) 


39.3 (28) 


47.7 (34) 


28.39 


28.22 


P12 


16.5 (12) 


33.2 (24) 


5.6 (4) 


28.31 


28.14 


ERS 


5.0 (6) 


5.0 (6) 


5.0 (6) 


27.20 


27.03 



a Y 09Sm for the ERS fields and Y 105w for the HUDF/P12/P34 
fields. 

Table 1. The total exposure time (in ksec) is listed for each filter, 
with the number of individual exposures given in parentheses. The 
final columns gives the 6 a and 7 a magnitude limits in the J-band 
- the 6 a is the limit of our catalog for candidate selection, and the 
luminosity function has been computed using y-drops brighter 
than t a. All magnitudes are on the AB system, measured in a 
O'.'6-diameter aperture with an aperture correction applied. 

SPARSAMPLE100, which non-destructively reads the array 
every 100 seconds. These repeated non-destructive reads of 
the infrared array allow gradient-fitting to obtain the count 
rate ("sampling up the ramp") and the flagging and rejec- 
tion of cosmic ray strikes. For the HUDF and flanking fields 
(Programme GO-11563) there were two exposures per orbit, 
with each MULTIACCUM comprising 16 reads for a total 
duration of 1403 sec per exposure. For the ERS images (Pro- 
gramme GO/DD-11359) there were 3 exposures per orbit, 
each with 9 or 10 reads and with a total time of 803-903 s 
per exposure. In Table[I]we list the exposure time (and num- 
ber of exposures) for each field and each spectral band. We 
note that for the V-band images of the HUDF, two visits (8 
exposures) were severely affected by image persistence, and 
as in Bunker et al. (2010) we exclude these from our data 
reduction. 



2.2 Data Reduction 

The IRAF.STSDAS pipeline calwf c3 was used to calculate 
the count rate and reject cosmic rays through gradient fit- 
ting, as well as subtracting the zeroth read and fiat-fielding. 
We used MULTIDRIZZLE (Koekemoer et al. 2002) to com- 
bine exposures taken through the same filter in each point- 
ing, taking account of the geometric distortions and map- 
ping on to an output pixel size of 0'.'06 from an original 
0'.'l3pix _1 . This was the same scale as we used in our anal- 
ysis of the Hubble Ultra Deep Field WFC 3 images (Bunker 
et al. 2010) and corresponds to a 2 x 2 block-averaging of 
the GOODSv2.0 ACS drizzled images. 

The final frames had units of electrons/sec, and 
we take the standard ACS zeropoints for the UDF im- 
ages. For WFC3, we use the recent zeropoints reported 
on http://www.stsci.edu/hst/wfc3/phot_zp_lbn during 
February 2010, where the F098M Y-band has an AB magni- 
tude zeropoint of 25.68 (such that a source of this brightness 
would have a count rate of 1 electron per second), and zero- 
points of 26.27, 26.25 & 25.96 for F105W, F125W & F160W. 
We note that the information in the image headers of the 
earlier images released in September 2009 is slightly differ- 
ent by 0.1-0.15 mag, with zeropoints of F105W Yzp = 26.16, 
Jzp = 26.10 & H Z p = 25.81 (as used in Bunker et al. 2010). 

The new WFC3 images of fields UDF-P12 and UDF- 
P34 were reduced using the latest version of calwf c3 (29 
October 2009 release), and we also re-reduced the UDF im- 



ages originally presented in Bunker et al. (2010) - the ear- 
lier paper had used a reduction using XDIMSUM and fits to 
the geometric distortion due to MULTIDRIZZLE not then 
being available for the newly-commissioned WFC3. For the 
UDF and flanking fields, we used a pixel fraction of 0.6 to 
recover some of the under-sampling. In each of these three 
fields we survey 4.18arcmin 2 in all exposures, with another 
0.67arcmin 2 surveyed at half the maximum depth in each 
field. For the ERS project, we reduced the data for all 10 
pointings using the same technique as described in Wilkins 
et al. (2010a, our analysis of the first 6 pointings). Here we 
used a pixel fraction of 1.0 in multidrizzle, as we only had 6 
exposures; a smaller pixel fraction would not fully populate 
the output repixellated grid with such a small number of 
exposures (this pixel fraction was also used for the i_~-band 
imaging of P12, as only 4 exposures had been taken at the 
time of writing - a small subset of the total still to be ob- 
served). We then mosaicked together the 10 ERS pointings 
in each filter, using inverse- variance weighting for the overlap 
regions, producing a field of fairly uniform depth covering 
37arcmin 2 , with a further 8arcmin 2 going less deep. 

In our final combined J-band image, we measure a 
FWHM of « O'.'l for point sources in the field. As most 
high-redshift galaxies are likely to be barely resolved (e.g., 
Bunker et al. 2004, Ferguson et al. 2004) we perform pho- 
tometry using fixed apertures of 0'/6 diameter, and intro- 
duce an aperture correction to account for the flux falling 
outside of the aperture. This correction was determined to 
be « 0.2 — 0.25 mag in WFC3 from photometry with larger 
apertures on bright but unsaturated point sources. We note 
that the ff-band images display significant Airy diffraction 
rings around point sources. For the ACS images, the better 
resolution and finer pixel sampling require a smaller aper- 
ture correction of w 0.1 mag. All the magnitudes reported in 
this paper have been corrected to approximate total magni- 
tudes (valid for compact sources), and we have also corrected 
for the small amount of foreground Galactic extinction to- 
ward these fields using the COBE/DIKBE & IRAS/ISSA 
dust maps of Schlegel, Finkbeiner & Davis (1998). The op- 
tical reddening is E(B — V) = 0.009, equivalent to extinc- 
tions of Assoip = 0.012, A 105w = 0.010, ^125™ = 0.008 & 
A lfi0w = 0.005. 

The geometric transformation and image re-gridding 
produces an output where the noise is highly correlated, 
hence measuring the standard deviation in blank areas of 
the final drizzled image will underestimate the noise. To as- 
certain the true significance of object detections, we deter- 
mine the real noise using several different techniques. As 
in Bunker et al. (2010) we also produced a crude combina- 
tion of the individual flat-fielded images using integer-pixel 
shifts. While this was not used for our science (as the sig- 
nificant geometric distortions were not accounted for, and it 
did not address the under-sampling of the PSF as "drizzle" 
does), this output frame had the advantage that the noise 
properties were preserved and adjacent pixels were uncorre- 
cted. We measured the standard deviation of the counts in 
blank areas of sky in this shift-and-add mosaic, and we ver- 
ified that the noise (normalized per unit time) decreased as 
the square root of the number of frames combined. The lim- 
iting magnitudes found using these uncorrelated "true-noise 
frames" are in good agreement with the STScI HST/WFC3 
Exposure Time Calculator (ETC) - Table [2] presents our 2 a 
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ACS/WFC3 2a Detection Limits (AB magnitudes) 



Field ID 


Center (J2000) 


Area (sq.arcmin) 


&435iu 


"606tu 


*775uj 


z 85(Mp 


^098m/105tu 


a Jl25w 


HlQ0 w 


HUDF 


03:32:38.4 -27:47:00 


4.2 


30.3 


30.7 


30.6 


30.0 


29.65 


29.70 


29.67 


P34 


3:33:05.3 -27:51:23 


4.2 




29.9 


29.5 


29.7 


29.45 


29.58 


29.41 


P12 


3:33:01.9 -27:41:10 


4.2 




29.9 


29.6 


29.6 


29.16 


29.50 


28.23 


ERS 


3:32:23.6 -27:42:50 


37.0 


29.1 


29.1 


28.5 


28.4 


28.02 


28.39 


28.10 



a Y c,$ m for the ERS field and Y 105m for the HUDF/P12/P34 fields. 

Table 2. Summary of observations. All magnitudes are on the AB system, and measured in a 0'.'6-diameter aperture with an aperture 
correction applied to correct to approximate total flux (for compact sources). 



limits in a O'.'6-diameter aperture, with the aperture correc- 
tion applied. We also produced a noise model based on the 
detector gain, readout noise and Poisson counts of the mea- 
sured background (including the instrument dark current), 
and verified that our sensitivity was well within 10 per cent 
of the expected noise. Finally, we measure the correlated 
noise (the standard deviation of the background counts) in 
the drizzled image mosaics which we use for our source de- 
tection and photometry, and use the relations in equation 
A13 of Casertano et al. (2000) to introduce a correction fac- 
tor which depends on the output pixel scale and the size of 
the "droplet" in the drizzling procedure ( "pixfrac" ) . We gen- 
erally found good agreement (at the 0.05 mag level) with our 
sensitivity measurements using the true-noise frames, except 
for the HUDF data where the corrected drizzle noise under- 
estimated the true noise by 0.1-0.2 mag, perhaps because of 
the large number of frames combined with small sub-pixel 
shifts. We adopted the sensitivity measurements from the 
true-noise frame, having checked that consistent results were 
produced by the ETC, the noise model, and the noise in the 
drizzle frame corrected for pixel correlations. Our measured 
noise in the HUDF is in good agreement with Bouwens et 
al. (2010a), but we note that McLure et al. (2010) appear to 
be ~ 0.3 mag less sensitive (although we note that their 5 a 
magnitude limit in a O'.'4-diameter aperture appears not to 
have been corrected to total magnitudes with an aperture 
correction, unlike in Bouwens et al. 2010a). 

The optical HST ACS imaging comes from the Hubble 
Ultra Deep Field (Beckwith et al. 2006), and we used the 
publicly-available v, i, z reductions of flanking field UDF- 
P12 provided by the UDF05 team (Oesch et al. 2007). We 
reduced the v,i,z ACS data for UDF-P34 from the HST 
archive, using MULT1DRIZZLE to combine a large subset 
of the data comprising blocks of data taken close in time with 
similar telescope roll angles, again using an output 0'.'03 pixel 
scale. These subsets of drizzled images were then registered 
and combined with IRAF. imcombine. Our combined images 
were 4.8ksec in v-band, 10.6 ksec in i-band and 26.8 ksec in 
z-band. All the ACS images were then block-averaged 2x2 
and registered with our drizzled WFC3 frames. 



2.3 Construction of Catalogues 

To perform the candidate selection we used the SExtractor 
photometry package (Bertin & Arnouts 1996), version 2.5.0. 
Since we are searching for F-drops (objects clearly detected 
in the WFC3 J-band but with minimal flux in the Y-band 
and ACS images), fixed circular apertures 0'.'6 in diameter 
were 'trained' in the J-image, and running SExtractor in 
dual-image mode those apertures were used to measure the 



flux in the same locations in the Y-band image. The same 
procedure was repeated between J-band image and all the 
other ACS and WFC3 images with different filters. For each 
waveband we used a weight image derived from the expo- 
sure map. To identify the objects we set the SExtractor pa- 
rameters to have a lower limit of 5 contiguous pixels above 
a threshold of 2a per pixel (data were drizzled to a scale 
of 0'.'06 pixel -1 ). We corrected the aperture magnitudes to 
approximate total magnitudes with the aperture correction 
appropriate for that filter. With this cut we were able to 
detect all significant sources, along with some spurious de- 
tections just above the noise limit or due to diffraction spikes 
from stars. We also impose a 6 a limit on the J-band mag- 
nitude for all fields, with the Jab magnitude limit listed in 
the last column of Table [I] Table [3] presents our photometry 
of K-drops from SExtractor, where we have corrected the 
magnitude errors returned by SExtractor for the effects of 
correlated noise in the drizzled images, using our "true noise 
frames" to determine the scaling factor (typically SExtractor 
underestimated the magnitude errors by a factor of fa 1.5 for 
pixfrac=0.6 used in most of our data, and a factor of « 2.6 
for pixfrac=1.0 as used in the ERS and the H-band of P12). 



3 CANDIDATE SELECTION 

Identification of candidates is achieved using the Lyman 
break technique (e.g. Steidel et al. 1996), where a large 
colour decrement is observed between filters either side of 
Lyman-a in the rest-frame of the galaxy. At z > 6, the flux 
decrement comes principally from the large integrated opti- 
cal depth of the intervening absorbers (the Lyman-a forest). 

At 2 « 8 — 9 the location of the Lyman-a break is 
redshifted to ~ 1.1pm - the WFC3 Y 105w / ogSm and J125™ 
are suitably located such that a 7.6 < z < 9.8 star forming 
galaxy will experience a significant flux decrement between 
these two filters (see Figure [T] & Figure [5j, although the 
selection efficiency drops at the extremes of this range. 

3.1 Contamination 

Searching for distant galaxies using only broadband pho- 
tometry means that contamination is a potentially serious 
issue. There are two main sources of contamination: objects 
whose intrinsic colours are similar to those of the target pop- 
ulation; and faint objects with intrinsically different colours 
but whose observed colours scatter into our selection be- 
cause of photometric noise. We note that the effect of tran- 
sient phenomena is not significant for the selection of Y- 
drops, since the WFC3 Y , J & H images were taken close 
in time. This is unlike our selection of z'-drops (e.g. Wilkins 
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3.6 3.8 4 4.2 

log 10 (X[A]) (observed frame) 



Figure 1. Top panel - Model (from the Starburst99, Leitherer 
ct al. 1999) spectral energy distribution (SED) of a redshifted 
2 = 8 star forming galaxy. Middle panel - Potential contaminants: 
Observed SED of a low-mass dwarf star (class: T4.5, Knapp et 
al. 2004) together with the model (Starburst99) SED of a 3.5Gyr 
Single-aged Stellar Population (SSP) at z = 2.5. The bottom 
two panels show the transmission functions of the combination of 
filters available to each field. 

et al. 2010a, 2010b; Bunker et al. 2010) where the ACS z - 
band and WFC3 Y-band were separated by many years, so a 
transient such as a supernova or high-proper-motion object 
which entered the Y-band but was absent at that location in 
the ACS could be erroneously identified as a Lyman-break 
galaxy. Indeed, a probable supernova was identified in the 
WFC3 imaging of the HUDF (e.g., Bunker et al. 2010). 

3.1.1 Intrinsically Red Objects 

There are two distinct types of objects whose apparent 
^i05W098m — J\25w colours are similar to those of Lyman- 
break galaxies at z as 8: lower-redshift (z ~ 2) galaxies 
have the Balmer/4000A break feature between the two fil- 
ters used, Y W5w /os$ m and J12&W , while some low mass dwarf 
stars, especially those of L and T spectral class, have low 
temperatures and broad absorption features that can mimic 
a spectral break. 

Examples of the spectral energy distributions (SEDs) 
of each of these types of object (a model 3.5 Gyr old single- 
aged stellar population at z — 2.5 and a T4.5 dwarf star) are 
shown in Figure [l] In the case of lower redshift galaxies the 
slope of the SED longward of the spectral break (i.e. long- 
ward of Y 10 5„/o98m) is somewhat redder than that predicted 
for a high-z star forming galaxy. The addition of a further 
filter at wavelengths redder than the J125™ filter (Hieow in 
this case) can then be used to discriminate between high-z 
and lower redshift galaxies (Figure |2J. L and T dwarfs con- 
tamination in the HUDF and P34 field is mostly ruled out 
by the Yiosm — J125™ colour selection we adopted. The ad- 



dition of Hiqow photometry is still important in excluding 
these objects in the ERS field (see Figure |3b, where the dif- 
ferent Y-band filter used provides less good discrimination 
using Y — J colour alone. 

In Figures [2] and [3] the positions of both the interlopers 
and the tracks expected for high-redshift star forming galax- 
ies are shown in the ( J125™ — Hi 6 q w ) - (Y 105w / Q98m — J125™) 
colour plane. With the exception of the lowest temperature 
T dwarfs where the Yogsm filter is employed (the ERS field), 
these interlopers form a distinct locus separate from z w 8—9 
star forming galaxies with UV spectral slope index /3 < 0.0 
(where fx = A' 3 is used as a model of the UV properties of 
star forming galaxies). 

Using this analysis it is possible to design a window 
in (Yios w /098m — Ji25w) - (Ji25w — Hieow) colour - colour 
space that selects mainly high-redshift star forming galaxies, 
while eliminating known contaminant populations. For the 
HUDF/P12/P34 fields (i.e. where we have Y 105w imaging) 
this YJH selection selection criteria is: 

(Yi05w — Jl25w) > 0.9 

(Yio5iu — Ji25«j) > 0.73 x (Ji25u> — Hieow) + 0.9 

(J1251D — -Hl60m) < 1.5 

The use of an alternative Y filter (Yiggm) in the ERS 
field necessitates the use of a slightly different criteria: 

(Y)98m — Jl25w) > 0.9 

(Yl98m — Jl25w) > 0.64 X (Jl25u> _ HlBOw) + 1-28 
{Jl25w — Hl60w) < 0.8 

We have designed our selection criteria to reject 
all known interlopers, while selecting most z w 8 — 
9 star- forming galaxies. Other groups have used similar 
colour:colour selection, but with slightly different colour 
cuts (e.g., Bouwens et al. 2010a). Although this may af- 
fect the surface density of candidates (due to a slightly dif- 
ferent redshift range and spectral range of spectral slopes 
probed for the LBGs, and a different contaminant fraction), 
the inferred luminosity functions should be similar as these 
selection effects are corrected for in the effective volume 
calculation. The window we obtain with such criteria ex- 
cludes a hypothetical population of z ~ 8 — 9 galaxies with 
Ji25jd — Hieow ^ 1 colours. Such a population would have 
extremely red UV spectral slopes, possibly due to massive 
dust reddening. The potential bias of our selection criteria 
and a more general analysis of the UV properties of the can- 
didates presented in this work is discussed in more detail in 
Wilkins et al. (2010c) where we conclude that the distribu- 
tion of UV spectral slope indices is consistent with being 
blue, with /3 ~ —2. 

3.1.2 Photometric Scatter 

Even with the selection criteria described above, we cannot 
prevent some objects being scattered into our selection win- 
dow because of photometric noise. At low signal to noise 
ratio, this contamination could be significant, and we im- 
pose another criterion to exclude those objects. To do that 
we use the deep optical imaging available in the ACS 6435^, 
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Figure 3. J125™ - Hieow and Y gg m - J\2hw colour - colour fig- 
ures for the ERS field showing our YJH colour selection window 
(grey shaded area), the location of our candidates, the predicted 
paths taken by high-rcdshift galaxies (solid lines, /3 = —3.0 , left, 
and /3 = 0.0 , right) and the location of possible contaminating 
sources. Contaminating sources include galactic stars (denoted by 
filled circles) and a passively evolving instantaneous burst of star 
formation that occured at 2 = 10 (dashed line). High-z candi- 
dates are denoted by black circles (where the size of the circle is 
an indication of the apparent Jab magnitude, with the brighter 
sources being bigger circles). Limits and errorbars are la. Ob- 
jects ERS.YD7 and ERS.YD8 are very marginal and hence not 
included in this figure (see Section 3.2). 
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Figure 2. J±25w - H 1(i0w and Y w5w - Jt25w colour - colour fig- 
ures for the HUDF (top) and flanking fields P34 (middle) and 
P12 (bottom), showing our YJH colour selection window (grey 
shaded area), the location of our candidates, the predicted paths 
taken by high-redshift galaxies (solid lines, = —3.0 , left, and 
jS = 0.0, right) and the location of possible contaminating sources. 
Contaminating sources include Galactic stars (O - T dwarf stars, 
with L and T stars being redder, denoted by filled circles) and 
a passively evolving 'early-type' galaxy (modelled as an instan- 
taneous burst of star formation at z = 10 followed by passive 
luminosity evolution, denoted by the dashed line). High-redshift 
candidates arc denoted by black circles (where the size of the cir- 
cle is an indication of the apparent Jab magnitude, with bigger 
circles indicating brighter sources). Limits and error bars are 1 a. 



«606uj, *775m and Zgsoip bands: because of the strong Lyman- 
a forest absorption, 2 « 8 — 9 galaxies should not have any 
significant flux in the 6435™, «606™ an d 1775™ bands, so we 
impose an additional bvi non- detection criteria for the se- 
lection of our candidates. All objects with a > 2<r detection 
in any of the b ii5w , «606™ and 1775^ are classified as con- 
taminants. The depths of these ACS images are given in 
Table 2. We note that the 2850™ filter does have a red tail 
which overlaps with the F 98m/i05iu-band filters used, so it 
is conceivable that a V-drop towards the lower end of the 
redshift selection might have residual 2-band flux. However, 
we found only one Y-drop (P34.YD5) with a ~ 2 a detection 
in 2-band. 

3.2 Candidate Galaxies at 2 m 8 — 9 

After imposing our selection criteria we are able to compile a 
list of candidate 2 sa 8— 9 star forming galaxies in the HUDF, 
UDF-P34, UDF-P12 and ERS fields. In Table|3]we list posi- 
tions and photometry of these objects, while thumbnails of 
the bvizYJH images of these candidates (where available) 
are presented in Figure [4] In total we find 24 Y-drop candi- 
dates (HUDF:6, UDF-P34:7, UDF-P12:2, ERS:9) covering 
a range of apparent Jab magnitudes of 27.0 — 28.5. In the 
three deep single WFC3 pointings, the number of candidates 
is fairly consistent from field to field, with 3, 4 and 2 Y-drops 
for the HUDF, UDF-P34 and UDF-P12 fields, respectively, 
at Jab < 28.2. 

There are 9 (of the 24) objects in the Y-drop list (Ta- 
ble^ which we flag as being more marginal than the other 
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candidates as they sit at the limits of our selection, although 
they are plausible z w 8 — 9 galaxies (our effective volume 
calculation already corrects for galaxies excluded as lying 
just outside the selection region). Candidates ERS.YD7 and 
ERS.YD8 in the ERS are flagged, as we only have a lower 
limit on the (Y — J) colour (they are < f a in Y"-band). 
Adopting the 1 a lower limit on the (Y — J) colour places 
them in or above the 'contaminant' triangular region of Fig- 
ure]^ fully consistent with entering our selection area. Simi- 
larly, objects ERS.YD2, ERS.YD5, ERS.YD9 and P34.YD7 
are nagged: using a f a lower limit on the (Y — J) colour 
these candidates would fully meet our selection criteria (see 
Figures [2] and [3), while a more conservative 2 a lower limit 
could potentially locate them just below our selection box, 
although with colours consistent with falling within the se- 
lection window. Deeper Y-band imaging is required to show 
unambiguously that they are not in the 'contaminant' region 
of the colourxolour space. Object P34.YD5 in P34 is also 
flagged, because it has a ~ 2 a detection in the z-band. There 
are no detections in «-, i- and Y-bands, though, so it is still 
a likely high-redshift (2 > 6) object - the z-band flux might 
be statistical fluctuation or perhaps a high-equivalent-width 
emission line within the z-band. 

We also flag as marginal two potential high redshift 
galaxies in field UDF-Pf2, on the grounds that the short ex- 
posure time of the //-band image in this field (TableJT]) made 
it impossible to measure the Hab magnitude. The upper 
limits on the (J — H) colours place them away from the red 
contaminant region with (J — H) > 1.5 (Figure [2J, but we 
require the UV luminosity in the //-filter (uncontaminated 
by the effects of Lyman-cv forest absorption) to infer the ab- 



solute UV magnitude (as described in Section 4.1 1. We now 
consider whether these single-band detections might be due 
to transients (such as was the case for the likely supernova 
in the WFC3 images of the HUDF, object zDO in Bunker et 
al. 2010). The P12 field was observed in J-band in two ob- 
serving blocks, with 8 frames taken on U.T. 2009 November 
02, and the other 16 frames taken over U.T. 2009 November 
10-15. As a check, we combined the two different epochs 
separately with "multidrizzle" . The magnitude of P12.YD1 
is consistent between the two epochs, with J = 28.07 ± 0.25 
(4.3 a) and J = 27.95 ±0.16 (6.8 a) respectively. However, 
P12.YD2 might show some variability in the J-band with 
J = 27.36 ± 0.13 (8.3 sigma) for the first block of data and 
J = 28.14±0.19 (5.8 sigma) for the second. Hence it is plau- 
sible that P12.YD2 might be a transient rather than a high- 
redshift Y-drop. When this WFC3 program (GO-11563) is 
complete, the //-band will be much deeper on P12, allowing 
a further check on the robustness of the candidates in this 
field. However, the two candidates in P12 represent less than 
10 per cent of our Y-drop sample, so will not quantitively 
affect our conclusions; for the moment, we exclude this field 
from our fitting of the UV luminosity function. 

3.3 Comparison with Other Studies 

We now compare our new list of candidates within the 
HUDF field with other groups' previous studies (Oesch et 
al. 2010, Bouwens et al. 2010a, McLure et al. 2010, Yan et 
al. 2010 and Finkelstein et al. 2010), and particularly with 
our previous paper (Bunker et al. 2010). A matched catalog 
between the Bunker et al. (2010), McLure et al. (2010) and 
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Figure 4. 2'.'4 X 2'.'4 (b)vizYJH thumbnail images of objects 
meeting our selection criteria in the analyzed fields. Within each 
field they are ordered by i/-band magnitude (brightest at the 
top). Objects marked with * are not more marginal candidates, 
see Section 3.2. The fields UDF-P12 and UDF-P34 do not have 
ACS 6-band imaging. 
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ID 


RA 


Dec 


YvB 


Jab (significance, a) 




(Y - J) AB 


(J - H) AB 


HUDF.YD2 


03:32:37.796 


-27:46:00.12 


> 29.76 


28.08 ± 0.12 


(8.9<t) 


28.20 ±0.14 


> 1.66 


-0.12 


HUDF.YD1 


03:32:42.873 


-27:46:34.58 


29.25 ±0.37 


28.15 ± 0.13 


(8.3cr) 


28.15 ± 0.13 


1.10 


0.00 


HUDF.YD3 


03:32:38.135 


-27:45:54.03 


29.79 ±0.61 


28.18 ± 0.13 


(8.1cr) 


28.10 ±0.13 


1.61 


0.08 


HUDF.YD4 


03:32:33.130 


-27:46:54.49 


29.85 ±0.65 


28.32 ± 0.15 


(7.1cr) 


29.29 ± 0.38 


1.53 


-0.97 


HUDFYD8 


03:32:43.082 


-27:46:27.75 


29.75 ±0.59 


28.45 ± 0.17 


(6.3<t) 


28.84 ± 0.25 


1.30 


-0.39 


HUDFYD9 


03:32:36.360 


-27:46:23.35 


29.78 ±0.61 


28.50 ±0.18 


(6.I0-) 


28.43 ±0.17 


1.28 


0.07 


P34YD1 


03:33:00.340 


-27:51:20.97 


29.71 ±0.68 


27.94 ± 0.12 


(9.1cr) 


28.27 ± 0.19 


1.77 


-0.33 


P34YD2 


03:33:09.657 


-27:51:16.45 


29.89 ±0.81 


27.95 ±0.12 


(9.0ct) 


27.61 ± 0.10 


1.94 


0.34 


P34YD3 


03:33:07.474 


-27:51:14.89 


29.85 ±0.78 


28.07 ± 0.14 


(8.I0-) 


28.86 ± 0.33 


1.78 


-0.79 


P34YD4 


03-33 04 857 


-27-51 -38 28 


on -jfi -|- n 


27.89 ± 0.11 


(9.5<t) 


27 83 -I- 1 3 


1.47 


0.06 


P34YD5* 


03:33:06.596 


-27:52:48.95 


> 30.20 


28.29 ± 0.17 


(6.6<t) 


28.20 ±0.18 


> 1.91 


0.09 


P34YD6 


03:33:06.901 


-27:51:32.54 


29.69 ±0.67 


28.35 ± 0.18 


(6.2a) 


28.03 ±0.15 


1.34 


0.32 


P34YD7* 


03:33:09.286 


-27:51:32.22 


> 29.91 


28.38 ± 0.18 


(6. Iff) 


28.85 ± 0.32 


> 1.53 


-0.47 


P12YD1* 


03:33:03.034 


-27:41:57.00 


29.64 ± 0.84 


28.01 ± 0.14 


(7.9cr) 


> 28.98 


1.63 


< -0.97 


P12YD2* 


03:33:00.545 


-27:41:46.80 


29.25 ±0.59 


28.03 ± 0.14 


(7.80-) 


> 28.44 


1.22 


< -0.41 


1 , 1 V.J . 1 1 J 1 


uo.oz .zo.ouy 


97-4Q-9fi ^Q 
-Z 1 .fto.ZU.OO 


ZO. / 1 


26.98 ±0.15 


(7.3<t) 


97 07 _i_ n 




n so 
u.oy 


njivD. i uz 


uo.oz.uz.you 


97-/1Q.^i nc: 
-z i .fto.oj- .yo 


\ 90 QQ 


26.98 ±0.15 


(7.3<t) 


9fi Qfi xn 17 


> 1 41 


12 


ERS.YD3 


03:32:29.790 


-27:43:01.09 


> 28.77 


27.03 ±0.16 


(7.0cr) 


27.83 ± 0.42 


> 1.74 


—0.80 


ERS.YD4 


03:32:40.904 


-27:40:12.37 


> 28.77 


27.06 ±0.16 


(6.80-) 


27.96 ± 0.47 


> 1.71 


-0.90 


ERS.YD5* 


03:32:18.414 


-27:43:45.99 


> 28.77 


27.06 ±0.16 


(6.8<r) 


27.35 ± 0.27 


> 1.71 


-0.29 


ERS.YD6 


03:32:05.022 


-27:45:53.93 


28.61 ± 0.91 


27.19 ±0.18 


(6.I0-) 


27.38 ±0.28 


1.42 


-0.19 


ERS.YD7* 


03:32:41.676 


-27:41:27.48 


> 28.03 


26.95 ±0.14 


(7.60-) 


26.56 ±0.13 


> 1.08 


0.39 


ERS.YD8* 


03:32:37.927 


-27:42:20.78 


> 28.14 


27.14 ± 0.17 


(6.3cr) 


26.89 ±0.18 


> 1.00 


0.25 


ERS.YD9* 


03:32:27.014 


-27:44:31.29 


> 28.57 


27.20 ±0.18 


(6.0ct) 


27.85 ± 0.43 


> 1.37 


-0.65 



Table 3. Y-band drop out candidate z 8 — 9 galaxies meeting our selection criteria in the JfUDF, P12, P34 and ERS fields. Objects 
are divided by field and then ordered by apparent Jab magnitude. Where the measured flux is < 1 a, an upper limit at the 1 a level is 
quoted. The significance of the J-band detection is shown in parentheses after the magnitude. Objects marked with * are more marginal 
candidates, see Section 3.2. The ERS field uses the narrower Y-band filter F098M, and the other fields use F105W. For the P34 field, 
we introduce a small offset to the astrometry in the HST image headers (given in Lorenzoni et al. http://arxiv.org/abs/1006.3545vl) of 
ARA = 0!03 to match the astrometry from GSC-2 and 2MASS. 



Bouwens et al. (2010a) samples has already been presented 
in Bunker & Wilkins (2009). 

Our refined HUDF sample, based on a new reduction 
of the HUDF data, has 6 Y"-band drop-outs. In Bunker et 
al. (2010) we presented a list of 7 Y-drop candidates within 
the HUDF field, the brightest four (in J-band) of which 
are reproduced with the new selection (HUDF-YD1,2,3 & 
4). Of the 3 other Y-drops from Bunker et al. (2010), one 
(YD5) has a discrepant (Y105™ — J125™) = 0.2 colour in the 
new data reduction, much bluer than our selection criteria 
of (Yio 5w — Ji25w) > 0.9. The faintest Y-drop in Bunker 
et al. (2010), YD7, is marginally too faint (J = 28.65) in 
our new reduction of the HUDF images to enter our new 
sample. However, applying our new colour selection criteria 
to the old photometry (where J = 28.44) would have re- 
sulted in the selection of YD7. The remaining one (YD6) 
is only marginally too blue for the Lyman-break selection 
in the newly-reduced data, with (Y/105™ — J/125™) = 0.89, 
very close to the (Y/105™ — J/i25u>) > 0.9 cut. This ob- 
ject has slight (~ 2 a) detections in the ACS bands, too, 
and does not meet the selection criterion (Yiosu, — Ji25iu) > 
0.73 x ( J125™ — Hi@o w ) + 0.9, so we did not include it in our 
list. Moreover, no other group has found or listed this object 
as a candidate. 

Two objects in our new catalog (HUDF.YD8 and 
HUDF.YD9) were not found in Bunker et al. (2010); our pre- 
vious study of Y-drops in the HUDF used slightly different 
magnitude and colour cuts (Jab < 28.5 and (Y — J)ab > 



1.0), and an older reduction and photometric zeropoints. 
These two objects were slightly too faint in the previous 
version of our HUDF reductions (J = 28.59 and J = 28.55, 
respectively) and slightly too blue ((Y/io5 m — J/125™) = 0.77, 
0.92 respectively) to be selected with our original criteria in 
Bunker et al. (2010). The new candidate HUDF.YD8 lies 
only 1 arcsec from the z-drop zD5 in Bunker et al. (2010), 
and it is conceivable that both objects might be physically 
associated and might have similar redshifts at z ~ 8. We 
note that no other group has identified HUDF.YD9 as a 
candidate. 

In Table [4] we show the Y-drop galaxy candidates from 
our HUDF catalog which have been previously reported with 
their corresponding catalog names from other groups, while 
in Table [5] we show all the objects found by these groups 
with colours or photometric redshifts compatible with being 
in our Y-drop redshift range, but which do not appear in our 
new catalog. We mark with a f the candidates that would 
be within our selection window if we adopt the photometry 
originally presented in the discovery papers, rather than re- 
measuring this with our new reduction of the HUDF WFC3 
imaging and the latest photometric zero-points. 

Most of the other HUDF candidates from different 
groups do not meet our selection criteria both because they 
are too faint in the J-band (class A in the Table [5| and 
because they are too blue, (Yiosu, — Ji25u>) < 0.9 (class B 
111 the Table |BJ. Only one candidate (z8-SB27 in Yan et al. 
2010) meets our selection criteria for brightness in the J- 
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BulO MIO YIO FIO Bunker '10 M10 F10 Y10 BolO Class 



ID BolO 



HUDF-YD2 


UDFy-37796000 


YD 2 


1939y 


zi 


S-B117 


200 


HUDF-YD1 


UDFy-42886345 


YD1 


1765y 


zi 


S-B092 


819 


HUDF-YD3 


UDFy-38135539 


YD3 


1721y 


zi 


S-B115 


125 


HUDF-YD4 




YD4 


2487 








HUDF-YD8 


UDFy-43086276 




2841y 


zi 


!-B088 


653 



Table 4. A list of Y-drops in the HUDF appearing in the cat- 
alogs of all previous analyses. We show in columns the different 
candidate ID used in this paper, in BolO (Bouwens et al. 2010a), 
BulO (Bunker et al. 2010), M10 (McLure et al. 2010), Y10 (Yan 
et al. 2010) and F10 (Finkelstein et al. 2010). 

band and the (Yiosm — J125™) colour, but is rejected on the 
basis of its location in the the J — H:Y — J colour:colour 
plane as a likely lower redshift Balmer-break galaxy (see 
Figure This galaxy is classified with letter 'C in the ta- 
ble. We note that Bouwens' candidate UDFy-37636015 (our 
YD7) has inconsistent photometry presented in Bouwens et 
al. (2009) and Bouwens et al. (2010a) - adopting the more 
recent photometric values from Bouwens et al. (2010a), YD7 
would enter our Y-drop selection (Table |5J. 

In summary, using our latest reduction of the WFC3 
images of the HUDF we are able to reproduce 4 of the 7 
Y-band dropout galaxies we first reported in Bunker et al. 
(2010); of two additional Y-drops in the new analysis, one 
has been reported elsewhere and one is a new discovery in 
the HUDF field. Remeasuring the photometry of Y-drop 
candidates presented elsewhere by other groups, we find that 
most would not enter our selection as they are too faint in 
J-band and/or are too blue in (Y — J), and hence are not 
as robust candidate z « 8 — 9 galaxies as our core sample. 



4 DISCUSSION 

4.1 The Luminosity Function of Y-drops at 

z^8-9 

From the observed surface density of Y-drops, as a func- 
tion of magnitude, we can recover the luminosity function 
of z ~ 8 — 9 galaxies in the rest-frame ultraviolet (observed 
by the WFC3 near-infrared niters). However, there is not 
uniform sensitivity over the redshift range probed by the Y- 
band drop-out technique; at the lower redshifts, the Y — J 
colour might be too blue to enter our selection, and at the 
higher redshift end of our range the effect of the Lyman-a? 
forest means that an increasingly large fraction of the J- 
band filter is absorbed, so only the most UV-luminous galax- 
ies will appear in our apparent-magnitude-limited sample. 
We quantify this effect, and hence constrain the luminosity 
function through our observed number counts. The probabil- 
ity of recovering a high-redshift galaxy as a function of red- 
shift and rest frame UV luminosity can be found with simu- 
lations. To perform these simulations we add into the images 
a large number of fake galaxies, with properties similar to 
those of the observed high-redshift population (i.e. compact 
with half-light radii thi ~ O'.'l, large Lyman-a forest decre- 
ment of Da ~ 0.99 and blue rest-frame UV colours). We 
then run our selection procedure and infer the probability of 
recovering such galaxies as a function of redshift and magni- 
tude (see Figure [5|. We adopt the effective volume approach 
as described in Steidel et al. (1999) and Stanway, Bunker & 



YD5t - - Z8-SD24 - AB 

YD6t - 

YD7t 2079y 213 z8-B114 y37636015 AB 

1107z" - - B 

1422 2055 z8-B041cf B 

800* - B 
3022 - AB 
640 - B 
z8-B094t 
z8-B087t 

z8-SB27 C 
z8-SB30 B 
z8-SD05 AB 
z8-SD02 AB 
z8-SD15 AB 
z8-SD52 AB 

f Object that would meet our selection criteria, assuming original 

photometry. 

* Object 800 appears in versions 1,2 & 3 of the arXiv:0912.1338| 
version of Finkelstein et al. (2010), but is absent from version 4 
and the Astrophysical Journal paper. 

** Object also found by Oesch et al. (2009b), named UDFz- 
44716442 and classified as a 2-drop. We included this object in 
our table since McLure et al. (2010) quote a photometric redshift 
of 7.60, on the edge of our selection range. 

Table 5. List of candidates from other studies of the HUDF (see 
Table |4j quoted by respective authors to be either a Y-drop or in 
our redshift range (7.6 < z < 9.8, see Section 3) but that we do 
not recover. Class A denotes objects too faint in J-band {Jab > 
28.51). Class B means that the colour selection criterion (Yiosuj — 
Jllhw) > 0.9. was not met. Class C is where the (J — H):(Y — J) 
colourxolour rules out selection (i.e., (Ylosto — J125™) > 0.73 X 
{Jl25w — H\qo w ) + 0.9 is not met). Objects marked f would meet 
our selection criteria if we adopt the original photometry (but 
which do not make our selection if we use the photometry from 
our new reduction of the imaging data) . 



McMahon (2003), such that the probability of recovering 
a galaxy in our survey depends on the redshift and abso- 
lute UV magnitudes, p(Muv, z), and from this the effective 
survey volume can be calculated (V e //). We use a Gaus- 
sian distribution of spectral slopes, with (/3) = —2.2 and 
cr(/3) = 0.5, reflecting the generally blue spectral slopes ob- 
served in Lyman-break galaxies at z > 6 (Stanway, McMa- 
hon & Bunker 2005; Bouwens et al. 2010b; Wilkins et al. 
2010c; Bunker et al. 2010). In Table [3] we have presented 
our list of candidate Y-drops, with colours consistent with 
being high redshift. These are good targets for spectroscopy, 
but in calculating the luminosity function we wish to restrict 
the sample to only the most reliable sources (to minimize 
biases through contamination by photometric scatter). In 
determining the luminosity function, we do not consider the 
P12 Y-drops, where the shallow H-band means we do not 
have secure if-band magnitudes. For the other fields, we in- 
clude only those galaxies from Table [3] detected at > 7 a in 
J-band which are not flagged as marginal; the only galaxies 
at > 7cr in J-band not included in the luminosity function 
fits are ERS.YD2 and ERS.YD7. 

We can then determine the best-fit luminosity function 
(Figure [bJ; we assume a Schechter (1976) functional form, 
where the number density of galaxies between L and L + SL 
is: 
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Figure 5. The probability of recovering simulated galaxies as a 
function of redshift and absolute rest-UV magnitude (Mi60o)- We 
have run simulations on all four of our fields. The mean redshift 
is denoted by a dot. 

<f)(L)dL = (j)*e~ x x° 1 

where x = L/L* . The Schechter function is parameterized 
by a faint end slope (a), a characteristic number density at 
the knee of the luminosity function (<j>*) and the characteris- 
tic luminosity, L* , corresponding to the absolute magnitude 
in the rest-frame UV (Mfjy, determined around 1600 A). 
Unfortunately we still do not have enough faint galaxies to 
constrain the faint end slope of this function, so we adopt 
three different values for the faint end slope, a = —1.5, 
— 1.7, —1.9, bracketting the value of a = —1.73 derived by 
Bouwens et al. (2006) for the i'-drops at z = 6 and for the 
z = 3 [/-drops (Reddy & Steidel 2009). We fit for the free 
parameters <f>* and M^y, and these are presented in Table[6] 

4.2 Evolution of the Luminosity Function with 
Redshift 

We now compare our measured best-fit luminosity func- 
tion parameters (Table 6) with previous estimates from the 
z w 8 — 9 galaxies in the HUDF alone. We also plot the 
z = 7 UV luminosity function from the z-drops of Wilkins 
et al. (2010b), and note that other LFs based on smaller 
datasets have derived similar parameters (e.g. McLure et al. 
2010; Oesch et al. 2010; Ouchi et al. 2010). Based on 5 Y- 
drops, Bouwens et al. (2010a) estimated M^ v — —19.45, 
assuming no evolution in <f>* and a from z ~ 6 (fixing <j>* — 
0.0011 Mpc -3 and a = -1.74). This is consistent with our 
determination of 0* = 0.00093 Mpc" 3 and M(j V = -19.5 
(where we have fixed a = —1.7 but fit both (f>* and M^y). 
We note that our measured charateristic number density is 
within « 20 per cent of the Bouwens et al. assumption, with 
M^jy nearly the same as the Bouwens et al. fit. McLure et 
al. (2010) suggest that the main luminosity function evolu- 
tion from z ~ 6 is in ci*, with My V broadly unchanged at 
Mjjy = —20, and <j>* z =§ ~ 5 x <j>l=a- However, this appears 
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Figure 6. The 2^8-9 rest-frame UV (1600A) luminosity 
function derived from HUDF, UDF-P34 and ERS WFC3 fields 
(circles, triangles and squares, respectively) together with con- 
temporary and lower-redshift comparisons. Solid lines denote the 
luminosity function at (2} = 3.05 (Reddy &; Steidel 2009) and 
(z) «= 7.0 (Wilkins et al. 2010b). 

to be marginally inconsistent with our number counts of Y- 
drops at the bright end - this parameter space is unavailable 
using the HUDF alone, but the larger volume we have in our 
current study enables us to fit both (j>* and M^y at z ~ 8 
- each of the 6 independent points in the luminosity func- 
tion (Figure |6| has w 3 galaxies in it (rather than 2 bins of 
2 — 3 galaxies previously). The McLure et al. (2010) pure- 
density-evolution scenario lies on the 1 a contour of our cj>* 
vs. L* reduced-^ 2 plot (Figure [71. Our results are entirely 
consistent with an evolution mainly in M^y since z — 7 
(and indeed since z = 3), with only a modest change in <f>* 
(consistent with no change in </}*). 

The integration of the luminosity function gives us the 
ultraviolet luminosity density, which is important for our 
purposes, since it is directly connected with the star forma- 
tion rate density. Figure [8] shows the UV luminosity density 
as a function of the magnitude down to which the luminos- 
ity function is integrated for each of the three best fits (one 
for each value of a adopted). 

4.3 The Star Formation Rate Density at z w 8 

We can use the observed H-band magnitudes of objects in 
the Y-drop sample to estimate their star formation rate 
from the rest-frame UV luminosity density around \ re st = 
1600 A. In the absence of dust obscuration, the relation be- 
tween the flux density in the rest-UV around this wave- 
length and the star formation rate (SFR in M0yr _1 ) is 
given by Ltjv = 8 x 10 27 SFRergss _1 Hz" 1 from Madau, 
Pozzetti & Dickinson (1998) for a Salpeter (1955) stellar 
initial mass function (IMF) with 0.1 M© < M, < 125 Af©. 
This is comparable to the relation derived from the models 
of Leitherer & Heckman (1995) and Kennicutt (1998). How- 
ever, if a Scalo (1986) IMF is used, the inferred star forma- 
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erg s 1 Mpc 
- 1 ) < 



-13 (> 0.01 M yr- 



^Mpc" 3 ]) 

x ) < -8 (> 10- 4 M Q yr- 



-19.34 
-19.5 
-19.66 



0.00117 
0.00093 
0.00070 



1.65 ( 0.0022 ) 
1.71 ( 0.0022 ) 
1.73 ( 0.0023 ) 



4.61 ( 0.0060 ) 
6.22 ( 0.0081 ) 
9.05 ( 0.0119 ) 



4.88 ( 0.0064 ) 
7.27 ( 0.0095 ) 
13.46 ( 0.0176 ) 



Table 6. The best fit values of M*^ 00 and <f>* f° r s Schechter function assuming fixed a £ { — 1.5,-1.7,-1.9} together with the UV 
luminosity densities (and star formation rate densities in parentheses) determined by integrating the luminosity function down to various 
limiting absolute magnitudes. 
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Figure 7. The significance contours for the reduced-x 2 fits of the 
Schechter UV luminosity function for the V-drops (solid contours 
on right, signifying reduced-x 2 = 1 an d 2 for the inner and outer). 
A faint end slope of a = —1.7 has been assumed, and the dot is the 
formal best fit, with M^ e00 = —19.5 (AB). The dashed contours 
(to the left of the z = 8 contours) denote the z = 7 luminosity 
function derived by Wilkins et al. (2010b) from the z-drops in the 
same WFC 3 fields as analysed here. The cross (on the left) is the 
z = 3 luminosity function for Lyman-break galaxies (Reddy & 
Steidel 2009). Evolution predominantly in M* is most consistent 
with the observational data. 



tion rates will be a factor of ~ 2.5 higher for a similar mass 
range. In the absence of a spectroscopic redshift, we assume 
that these lie at the predicted average redshift for galax- 
ies obeying our colour cuts (Figure |5J. For the luminosity 
functions considered, the predicted mean redshift is around 
(z) — 8.6 for a spectral slope j3 ~ —2 and Muv = —20.5. 
The i/-band probes the rest-UV above Lyman-a, and is 
unaffected by the forest absorption. Our 6 a limit for the 
HUDF is Hab = 28.5, equivalent to an absolute magnitude 

of M " — —18.6 at z — 8.6, corresponding to an inferred 
1600A ' ^ to 

star formation of 1.5 Mq yr _1 . This is equivalent to 0.1 -L*_ 3 
for M uv = -20.97 at z = 3 (Reddy & Steidel 2009) . 

With our measured Schechter luminosity function pa- 
rameters (0* = 0.00093 Mpc" 3 , M* = -19.5 assuming 
a = — 1.7 as at 2 = 3 and z = 6), the total star forma- 
tion rate density is 0.0022 Mq yr" 1 Mpc~ 3 integrating down 
to our luminosity limit of M uv ■ w -18.5 mag (AB). This 
should be regarded as a robust lower limit on the star for- 
mation rate density, as dust obscuration may affect the rest- 
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Figure 8. The UV luminosity density (left axis) and star forma- 
tion rate density (right axis) as a function of the rest-UV (M1600) 
absolute magnitude down to which the luminosity function is in- 
tegrated. We show our best fit luminosity functions at z ~ 8 — 9 
assuming a = —1.5, —1.7, —1.9. The shaded grey box denotes the 
observed region with the remainder inferred from extrapolation 
of the luminosity function. The horizontal lines show the UV lu- 
minosity density required to reionize the Universe at this redshfit, 
assuning a clumping factor of C = 5, and an escape fraction of 
/csc = 0.5 (top line), 0.7 (middle) and 1 (bottom). 



frame UV continuum, and also galaxies fainter than our se- 
lection limit will also contribute to the integrated UV light 

density. If instead we integrate down to M uv ° = — 13 (cor- 

1600A v 

responding to O.OlM yr _1 ) and the total star formation 
rate density is 0.0081 Mq yr _1 Mpc -3 . These star formation 
rate densities are a factor of ~ 10 lower than at z ~ 3 — 4, 
and even a factor of w 3 — 5 below that at 2 w 6 (Bunker et 
al. 2004; Bouwens et al. 2006). 

4.4 Implications for Reionization 

The ionizing UV photons produced by the most massive 
(OB) stars might be critical in reionization and keeping the 
Universe ionized at z ~ 6 — 11. However, work at z ~ 6 
has shown that under standard assumptions of the IMF, 
escape fraction and clumping of the gas, the observed pop- 
ulation of Lyman break galaxies produce insufficient flux 
down to AB ~ 28.5 mag (Bunker et al. 2004), and the "pho- 
ton drought" is even more severe at 2 w 7 (Wilkins et al. 
2010b). We now compare our measured UV luminosity den- 
sity at 2 w 8 — 9 (quoted above as a corresponding star 
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formation rate density) with that required to ionize the Uni- 
verse at this redshift. Madau, Haardt & Rees (1999) give the 
density of star formation required for reionization (assuming 
the same Salpeter IMF as used in this paper): 



PSFR : 



0.012 M yr _1 Mpc" 

/esc 



1 + z 
1 + 8.6 



fib hyp 
0.0462 



We have updated equation 27 of Madau, Haardt & Rees 
(1999) for a more recent concordance cosmology estimate 
of the baryon density from Larson et al. (2010), fib /iioo = 
0.022622. The reionization requirement at z ~ 8.6 is a factor 
of 2.5 times higher than that at z « 6, as the number of 
photons needed rises as (l + z) 3 . 

In the above equation, C is the clumping factor of neu- 
tral hydrogen, C = (phi) (phi)~ 2 - Early simulations sug- 
gested C w 30 (Gnedin & Ostriker 1997), but more recent 
work including the effects of reheating implies a lower con- 
centration factor of C ~ 5 (Pawlik et al. 2009). The escape 
fraction of ionizing photons (/ esc ) for high-redshift galaxies 
is highly uncertain (e.g., Steidel, Pettini & Adelberger 2001, 
Shapley et al. 2006), and it is possible that the escape frac- 
tion of ionising photons may be linked to the escape fraction 
of Lyman-Q photons (Stark et al. 2010), which may mean 
that high escape fractions could be tested through future 
line emission line searches with spectroscopy and narrow- 
band imaging. Even if we take the upper limit of / eS c = 1 
(no absorption by H i) and a very low clumping factor, the 
required total star formation rate density for reionization is 
0.012 M yr~ 1 Mpc~ 3 . This is a factor of ~ 5 higher than 
our measured star formation density at z w 8 — 9 from 
Y-drop galaxies brighter than Muv = —18.5 (our approx- 
imate limit). As shown in Table [6] and Figure |8j the re- 
quired UV luminosity density can only just be achieved (if 
/esc = 1) by integrating down to Muv = —13 (i.e., ex- 
trapolating the Schechter function to w 100 times fainter 
than our observed limit) and then only for a steeper faint 
end slope of a = —1.9 rather than a = —1.7. Adopting a 
less unrealistic value of / CS c = 0.7 (which is still high com- 
pared with observed values at lower redshift) the required 
total star formation rate density for reionization would be 
0.017 Mq yr _1 Mpc -3 , then the F-drop population can only 
provide sufficient ionizing photons if the faint end slope is 
very steep (a < —1.9) and the Schechter function is inte- 
grated down below Muv = —8 (corresponding to a star for- 
mation rate of only 10~ 4 Mq yr" 1 ). We note that recent the- 
oretical papers indicate that the reionization process itself 
may have been "photon-starved" (e.g., Bolton & Haehnelt 
2007), consistent with the extrapolation of our observational 
constraints. 

However, the assumption of a solar metallicity Salpeter 
IMF may be flawed: the colours of z ~ 6 i'-band drop-outs 
are very blue (Stanway, McMahon & Bunker 2005), with /3 < 
—2, and the recent WFC3 J- and ff-band images show that 
the z ~ 7 z'-drops also have blue colours on average (Bunker 
et al. 2010; Bouwens et al. 2010b; Wilkins et al. 2010c). 
Continuous star formation with a Salpeter IMF produces a 
UV spectral slope of ft m —2 if there is no dust reddening. 
The fact that we observe even more blue slopes than this 
(/3 < —2) could be explained through low metallicity, or a 
top-heavy IMF, which can produce between 3 and 10 times 
as many ionizing photons for the same 1600 A UV luminosity 
(Schaerer 2003 - see also Stiavelli, Fall & Panagia 2004). 



Alternatively, we may be seeing galaxies at the onset of star 
formation, or with a rising star formation rate (Verma et 
al. 2007), which would also lead us to underestimate the 
true star formation rate from the rest-UV luminosity. We 
explore the implications of the blue UV spectral slopes in 
z > 6 galaxies in a forthcoming paper (Wilkins et al. 2010c). 



5 CONCLUSIONS 

In this paper we have presented a search for galaxies at 
7.5 < 2 < 10 using the latest HSTWFC3 near-infrared data, 
based on the Lyman-break technique. Searching for galaxies 
which have large (Y — J) colours (Y-drops) on account of 
the Lyman-alpha forest absorption, and with (J—H) colours 
inconsistent with being low-redshift contaminants, we iden- 
tify « 20 candidates at redshift z w 8 — 9 over an area 
of w 50 square arcminutes. Our deepest field (the HUDF, 
covering 4.2arcmin 2 ) reaches Jab = 28.5 at 6 a, while the 
wide-area ERS data (comprising 10 WFC3 pointings cover- 
ing 37arcmin 2 ) reaches Jab = 27.2. The surface densities 
of candidates as a function of limiting magnitude appear 
broadly consistent between our 4 fields, although these all lie 
within 10 arcmin. Previous searches for Y-drops with WFC3 
have focussed only on the HUDF, and our larger survey has 
trebelled the number of robust Y-drop candidates, as well as 
providing a number of brighter Y-drops (with Jab ~ 27.0 
rather than Jab > 28.0 as in the HUDF). These brighter 
sources may be more amenable to spectroscopic follow-up. 

For the first time, we have a sufficient number of z « 
8 — 9 galaxies to fit <j>* and M* assuming a Schechter lumi- 
nosity function (previous estimates had to fix one of these 
parameters). We confirm that there is large evolution from 
z = 3, particularly in the bright end of the luminosity func- 
tion, in the sense that there are far fewer UV-bright galaxies 
at z w 8 — 9 than in the more recent past. There is also ev- 
idence for evolution from z — 6 — 7 to 2 = 8 — 9, with this 
being consistent with most of the change occurring in M* 
rather than </>*, with M* being fainter at higher redshift. We 
are unable to obtain a good constraint on the faint-end slope, 
a, which will potentially require deeper data over a wider 
field (as might be provided by NIRCAM on the James Webb 
Space Telescope). The candidate 2 w 8 — 9 galaxies we detect 
have insufficient ionizing flux to reionize the Universe, and 
it is probable that galaxies below our detection limit provide 
a significant UV contribution. However, adopting a similiar 
faint-end slope to that determined at 2 = 3 — 6 (a = —1.1) 
and a Salpeter IMF, then the ionizing photon budget still 
falls short if / osc < 0.5, even integrating down to Muv = —8. 
A steeper faint end slope and a low-metallicity population 
(or a top-heavy IMF) might still provide sufficient photons 
for star-forming galaxies to reionize the Universe, but con- 
firmation of this might have to await the James Webb Space 
Telescope. 
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